State tomography via weak measurements 
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Recent work has revealed that the wave function of a pure state can be measured directly and 
that complementary knowledge of a quantum system can be obtained simultaneously by weak mea- 
surements. However, the original scheme applies only to pure states, and it is not efficient because 
most of the data are discarded by post-selection. Here, we propose tomography schemes for pure 
states and for mixed states via weak measurements, and our schemes are more efficient because we 
do not discard any data. Furthermore, we demonstrate that any matrix element of a general state 
can be directly read from an appropriate weak measurement. The density matrix (with all of its 
elements) represents all that is directly accessible from a general measurement. 



For a group of blindfolded observers, who can each 
only touch one part of the elephant, to construct an accu- 
rate representation of an elephant, the group must com- 
bine information about different parts of the elephant 
(FIG la). In quantum mechanics, to have a complete 
description of a quantum system, in particular, a com- 
plete knowledge of a quantum state, we must combine 
information about complementary aspects (FIG lb). 

An ideal (strong) quantum measurement of a certain 
observable only gives the probabilities of obtaining the 
eigenvalues of the observable, and the statistical results 
of the measurements reflect the diagonal terms of the 
density matrix in the eigenbasis of the observable. Mea- 
surements of different (complementary) observables pro- 
vide the diagonal terms of the density state for different 
bases. If we perform measurements for a sufficient num- 
ber of bases, we can reconstruct the density state, which 
is the idea of state tomography. 

An important difference exists between the tomogra- 
phy of a quantum state and its classical counterpart, e.g., 
the description of an elephant by blindfolded observers. 
In this classical example, the description (of the elephant) 
by each man is not exclusive, and the different aspects of 
perception can be simply combined to produce an overall 
description of the elephant (FIG la). However, quantum 
physics forbids simultaneous knowledge of complemen- 
tary observables because precise knowledge of a certain 
aspect necessarily implies uncertainty for the complemen- 
tary aspects [l[. Therefore, one cannot simultaneously 
perform ideal measurements of complementary observ- 
ables in quantum physics (FIG 2a). To determine the 
quantum state of a system, different measurement se- 
tups, each for a particular observable, are required. To 
determine the general state of a d-dimensional quantum 
system, at least d + 1 different experimental setups are 
needed. Each setup must be devoted to the measurement 
of one of the complementary observables. 

An alternative state tomography approach is possible. 
Instead of obtaining maximum information of a particu- 
lar observable by an ideal measurement, we can perform a 
weak measurement Q • The state of the quantum system 
is only slightly changed during its weak interaction with 
the measuring device; therefore, weak measurements of 




FIG. 1. (a) Compatible aspects are combined in classical 
physics, (b) while complementary aspects are combined in 
quantum physics. 
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FIG. 2. (a) A strong measurement only reveals a certain 
aspect of the quantum state, (b) while a weak measurement 
enables us to perceive (measure) incompatible aspects of the 
quantum state simultaneously. 



a set of complementary observables can be performed 
simultaneously Q (FIG 2b). Thus, fewer experimental 
setups are needed for state tomography via weak mea- 
surements m. 

Weak measurements have been well described in the 
literature. Since they were first introduced by Aharonov, 
Albert and Vaidman weak measurements have been 
realised in experiments [5|-l7j], and have provided new in- 
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sights into the study of paradoxes and fundamental prob- 
lems in quantum theory [8j-[l2j. Weak measurement has 
also been used as a practical tool for amplifying weak 
signals and signal-to- noise ratios [I3I - I20I ]. Signal ampli- 
fication via weak measurement usually involves a pre- 
selection and a post-selection that are nearly orthogonal, 
and the original formalism of weak measurement [3, Hl[ 
is not sufficient for the phenomena in this regime because 
it only retains the first-order terms of the interaction 
strength [22Tj2^ ]. Extensions to the case of general pre- 
selection are given in 

El 

[26| . A framework that retains 
the high-order terms for general pre-selection and post- 
selection is given in [27j , and specific cases are further 
studied in [28l430| . A nonperturbative theory of weak 
pre- and post-selected measurements is given in (3lj | . The 
general results for the weak measurement of a pair of 
complementary observables are provided in [20| . Related 
concepts, such as the contextual values [32| and modular 
values [33j], were also introduced. Weak measurements 
can also be used for state tomography. Lundeen et al. [H 
directly measured the wave function of a pure state by 
a weak measurement followed by a post-selection. The 
scheme was extended to the mixed state by sequential 
weak measurements of pairs or triple products of com- 
plementary observables 34]. The results of the weak 
measurements could also be interpreted in terms of the 
Kirkwood-Dirac representation via complex joint prob- 
abilities (3f| H|| . Reviews and additional references on 
weak measurements can be found in [3l|, H3, H|[ ■ 

The results in [H demonstrate that the wave function 
of a pure state can be directly obtained from a single 
experimental setup, and that one can directly measure 
both the absolute values and the phases of the coeffi- 
cients of a pure state in a certain basis. However, the 
scheme proposed in ji[ applies only to pure states, and 
it is not efficient because most of the data are discarded 
due to the post-selection. Here, we propose tomography 
schemes for both pure states and mixed states via weak 
measurements. Our schemes are more efficient because 
we do not discard any data. We also show that any (di- 
agonal or off-diagonal) element of the density state can 
be directly determined as the average pointer shift in an 
appropriate weak measurement. 

Results 

The idea. In the scheme proposed in Q, the wave 
function to be measured is a continuous function of the 
position. The essence of the scheme is reviewed in finite 
dimensional Hilbert space as follows [13, HI] . Suppose a 
quantum system with a d-dimensional Hilbert space is in 
an unknown state 
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where the coefficients {ipi} in a certain basis {]£!,)} need 
to be determined. The scheme for measuring the coef- 
ficients {ipi} consists of a series of weak measurements 



of the observables Ai — |a,) {cti\, each followed by the 
same postselection, i.e., a projection onto the final state 
\bo) = -^2^2^=0 l a *)- The we& k value of the observable 
A4 is defined as 
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Both the real and the imaginary parts of the weak value 
Wi have physical meanings and can be determined exper- 
imentally because they correspond to the average shifts 
of the pointer position and the momentum, respectively. 
From Eq. ([2]), we know that the weak values Wi are di- 
rectly proportional to the coefficients ipi that we want to 



measure and that the factor 



can be determined 



by the normalisation condition (up to an unphysical over- 
all phase of \ip}). Therefore, we obtain a direct measure- 
ment of the coefficients ipi, thus a direct measurement of 
the wave function l^). 

The essential point in this scheme is the choice of the 
post-selected state |&o)- Its overlap with each basis state 
\di) has the same magnitude and phase; therefore, the 
factor K = tpi/Wi = Vd{bo\ip) does not depend on i 
and can be determined by normalisation. Because the 
change of the system state due to a weak interaction with 
the pointer is negligible, the success probability of post- 
selection is given by P = | (bo\tp) \ 2 ■ Therefore, only a 
fraction P of the data is retained, and the remainder is 
discarded due to the post-selection. When the dimension 
d is large, as in the case of a continuous wave function, 
the majority of the data is discarded. 

State tomography of a pure state. The part of the 
data that corresponds to the failure of the post-selection 
can be retained. We replace the final post-selection by 
a complete projective measurement onto the basis states 
{|6j)}. The inner products /3ji = (bj\a,i) are fixed when 
the two sets of the basis states are chosen. We organise 
the data according to the final state: if the final state is 



\bj), then the weak value of Ai = |aj) (a 
w = (bj\ai) (oj|V>) = (bj\ai) ^ = 



is given by 
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Therefore, for each fixed j, the weak values of different 
Ai give the relative ratios of the coefficients ipi'- 



(J) _ 



(bj\ ai 



(4) 



Because (bj\ip) does not depend on i, the coefficients ipi 
are directly proportional to -j 21 , where the weak values 
Wji are read-outs directly from the experiments and the 
constants /3ji = (bj\a,i) are fixed when the two sets of ba- 
sis states are chosen. The super index (j) in ipf"^ denotes 

the coefficients ipi that we obtain from the jth subset 

(7) 

of the data. Ideally, ip\ should not depend on j; how- 
ever, they may depend on j due to imperfections and 



3 



noise in the experiments. The magnitude of the factor 



(b 3 



can be removed by normalisation: 



(5) 



The phase factor tpj has no physical meaning because it 
corresponds to the overall phase of the pure state we are 
attempting to measure. This factor can be removed by a 
convention. For example, we require ipo > (or the co- 
efficient with the largest magnitude be positive). Using 
this approach, we retain all of the data and separate the 
data into d subsets according to the final states 
For the jth subset of the data that corresponds to the 
fixed final states \bj), we have an estimation of the coeffi- 
cients ipl^ from Eq. (0J or §5§ . A different subset of the 
data gives a different estimation for the set of coefficients 
{ipi}. The differences between the estimations indicate 
the amount of error or noise in the experiment. 

When the basis {|&j}} is mutually unbiased to the basis 

{|aj)}, i.e., /3ji — (bj\ai) — ^7=-, © can be simplified to 



s'We-^'Wtf/./El^'l 



(6) 



The original scheme proposed in [J] corresponds to the 
case in which \b ) = -7= J3i=o l a *)' anc ^ on ^ ^he data 
corresponding to the successful post-selection of |&o) are 
retained. Here, we retain all of the data, and each subset 
of the weak values corresponding to a fixed final state 
gives an estimation of the state to be measured. The 
probability of obtaining the final state \bj) is given by 
I (bj\ip) 1 2 because the change of the system state during 
the weak interaction is negligible. When the state to 
be measured is almost orthogonal to a certain final state 
\bj), then the relative frequency to obtain the final state 
\bj) in the post-selection is small. However, the corre- 
sponding weak value has a large magnitude, which could 
result in a large signal-to-noise ratio when the dominant 
noise is due to systematic error or imperfections. 

In the above scheme, weak measurements for a com- 
plete set of projective operators {Ai = |a») (di\} are 
needed. For the tomography of a pure-state wave func- 
tion \ip), weak measurements of a single observable in- 
stead of a set of observables are sufficient. This single 
observable can be chosen as follows (an alternative choice 
is presented in the Methods). We consider a pure state 
\tp) such that its overlap with each postselected state is 
nonzero, i.e., (bj\(p) for each j. We perform a weak 
measurement of P v = \(p) (<p\. When the post-selected 
state is \bj), the weak value of P v is defined as 



Wj = 



(bjlPM _ (bAtp) (^#) 



Therefore, 



(bj\<P) 



(¥#) = Vj 



(7) 



(8) 



Here, r\ 3 — is completely determined from the ex- 

perimental data and our choice because Wj is directly 
obtained from the experiment and both \<p) and 
are fixed by our choice. From Eq. ([5]). we can recon- 
struct the state \ip): 



(9) 



where K is determined by the normalisation condition 
up to a phase. 

State tomography of a mixed state. To this point, 
we have only considered tomography schemes for a pure 
state. If the state we like to measure is a general mixed 
state p, we also have a tomography scheme via weak mea- 
surements. For this case, we perform weak measurements 
for a complete set of projective operators Ai = \at) (aj|, 
followed by a final projective measurement onto the basis 
For the initial state p and the final post-selected 
state \bj), the average shifts of the pointer position and 
the momentum in a weak measurement of A4 correspond 
to the real and imaginary parts of the weak value [23, W\ 
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H\bi)(bj\p) 



(bj\p\bj 



From Eq. (|10p , one can express the matrix elements of p 
either by the basis {|aj)} 



\p\ a j) = ^2ih\p\h 



(6]fe|ai) ^ Pki 



or by the basis {!&■,■)} 
(b i \p\b j ) = (b j \p\b j )J2w jk 



(11) 



(h\cik) 
(bj\a k ) 



( 12 ) 

where the probability Pj = (bj\ p \bj) corresponds to the 
frequency of obtaining the final state \bj) in the experi- 
ment. The weak values Wji and the probabilities Pj are 
directly accessible from the experiment. The relation- 
ships in Eqs. pT|) and (TT^|) show that the matrix elements 
of the mixed state can be written as linear summations 
of the weak values that are directly accessible from the 
experiment. 

Partial state tomography. Using the idea of weak 
measurements, we can also selectively measure some of 
the matrix elements (diagonal and off-diagonal terms) of 
a state directly. This is especially useful when we are 
only interested in one or a few matrix elements of the 
state and we do not need to perform state tomography 
for the entire density matrix. Suppose we are interested 
in measuring a particular matrix element (a\p\b). We 
consider two different cases according to whether | a) and 
I b) are orthogonal. For the first case, suppose (b\a) 7^ 
0. Here, we can simply perform a weak measurement of 
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| a) (o|, with a post-selection onto the state \b). The weak 
value is defined as 



W = 



(b\a) (a\p\b) 
(b\ P \b) ' 



Thus, 



(a\p\b) 



(b\ P \b) 
(b\a) 



W. 



(13) 



(14) 



W is the weak value, and (b\p\b) is the probability of 
success for post-selection. Both values are directly ob- 
tained from the experiment. (b\a) is only a fixed factor. 
Therefore, (a\p\b) is determined by Eq. fj 14|) . For the 
second case, suppose that (b\a) = 0. We then perform 
a weak measurement of the observable C = |c) (c| with 
|c) = 7^(| a ) + \b)) and a follow-up projective measure- 
ment onto a basis including both \a) and \b) as the basis 
states. We must only retain the data when the final state 
is either \a) or \b). If the final state is \a), then the weak 
value of C is given by 



= (a\c) (c\p\a) = 1 
(a|p|a) 2 { 



(b\p\ 



(a\p\a) 



)• 



(15) 



If the final state is \b), then the weak value of C is given 
by 



(b\c) (c\p\b) = 1 

(6|p|6> ~ 2 l + (6| P |6) J - 



We have 



(6|p|o)=<a|p|a) (2W-1) 
(a\p\b)={b\p\b) (2W-1) 



(16) 



(17) 
(18) 



The post-selection probabilities {a\ p \ a) and (&| p \ b) and 
the weak values W and W are directly accessible from 
the experiment; therefore, we can determine the matrix 
element (a\p\b) from either Eq. (jTTJ) or (|18l) . which 
should agree with each other. Any discrepancies between 
them could be used as indicators of the noise and imper- 
fections in the experiment. 

Experimental scheme. Here, we discuss how our 
schemes can be realised experimentally. For the weak 
measurement of an observable A, the interaction between 
the measuring device and the system is generally mod- 
elled by a Hamiltonian T~L — gA eg) p8(t — to), where p de- 
notes the momentum of the pointer (measuring device). 
Because the interaction in a weak measurement does not 
change the state of the system significantly, weak mea- 
surements of several observables (whether they are com- 
mutative or not) can be performed simultaneously. The 
interaction of the simultaneous weak measurements of a 
set of observables {A{\ can be introduced by the Hamil- 
tonian: 



n = ^2g l A l <E) Pt 5{t-t ). 



(19) 



Here, pi denotes the momentum of the ith pointer that 
is coupled to the observable A^. {Ai} is an arbitrary 
set of observables of the system. They could be a set 
of commuting observables (for example, Ai = \a,i) (a,i\) 
or even a set of complementary observables. Suppose 
the initial state of the system is p s , and the initial state 
of the pointers is pd = pi ® pi ® • • • , where pt denotes 
the initial state of the ith pointer. Suppose the post- 
selection of the system is a projection H, (for example, 
Hj = \bj) (bj\). Then the average position shift Sqi and 
the average momentum shift dpi of the ith pointer are 
given by (see the derivation in the Methods) 



8qi = gtReWji 

Spi = 2g l ImW jl (Ap i ) 2 . 



(20) 
(21) 



The weak value Wji of Ai, for a general initial state p s 
and a general postselection H,, is given by [2(| 



Wji = 



Tr(H 3 A lPs ) 
Tr(U jPs ) ■ 



(22) 



The shifts of different pointers can be read out simulta- 
neously, while either the position shift or the momentum 
shift is recorded at each time for each pointer. With the 
same Hamiltonian, we can obtain all of the weak values 
Wji. From these weak values, we can obtain all of the 
elements of the density matrix p s according to the state 
tomography strategy proposed in this article. 



Example of application. The state tomography 
schemes proposed here can also be used to detect a tiny 
parameter encoded in a quantum state. For example, to 
detect a tiny phase difference •& picked up by two orthog- 
onal states of a qubit (a photon) that passes through a 
certain medium (a birefringent crystal), we could pre- 
pare a qubit in the initial state ^= (|0) + |1)), and seek 
to determine the tiny phase d encoded in the outgoing 
state A= (|0) + e l ® |1)J. To determine we can per- 
form a weak measurement of the observable |1) (1| by 
introducing an interaction H = g |1) (1| ® pS(t — to) be- 
tween the qubit system and a measuring device and a 
follow-up post-selection of the qubit system onto the 
state ^ (|0) - |1)). The weak value is W = 1 - 
and the average shift of the pointer momentum is given 
by Sp = — 2g(Ap) 2 i. Although the phase t? is small, 
the average shift of the pointer momentum is not nec- 
essarily small, and this has certain advantages over the 
quantum-entanglement based scheme in Ramsey inter- 
ferometry [ioj|, which relies on the use of GHZ states of 
many qubits, which are an expensive resource. 

Discussion 

In this article, we have proposed several schemes for state 
tomography via weak measurements. An advantage of 
our schemes is that we need to measure fewer observ- 
ables compared with the standard method of state to- 
mography based on ideal measurements. For example, 



5 



a simple standard state tomography requires ideal mea- 
surements in at least d + 1 different bases; therefore, one 
needs to use at least d+1 different experimental setups. 
Here, our tomography strategy requires weak interaction 
in one basis ({|ai)}) and post-selection in another basis 
and the same interaction Hamiltonian is used 
for all data collection. 

Compared with the scheme in our strategy ap- 
plies not only to pure states but also to mixed states. 
Our strategy is more efficient because we retain all of 
the data and thus reduces the number of experimen- 
tal runs needed. Compared with the schemes in [34{ in 
which a density matrix is directly measured via sequen- 
tial weak measurements of pairs or triple products of ob- 
servables, our schemes are based on a single-time weak 
measurement followed by a strong complete projective 
measurement. Our interaction Hamiltonian in Eq. (|19[) 
is straightforward and differs significantly from those in 
the schemes of [34| • In the previous schemes d, [34| , the 
observables A4 = |etj) (aj| require separate weak measure- 
ments at each time for a particular i. In our schemes, the 
set of observables {Ai} can be measured simultaneously 
via the interaction introduced by the simple Hamiltonian 
in Eq. (fTO]). 

Interestingly, any matrix element of a density state is 
directly accessible from a suitable weak measurement. 
Ideal measurements give probabilities and expectation 
values, which are only linear combinations of the ma- 
trix elements of the density operator. Based on all of 
these facts, it is natural to draw the conclusion: what is 
and can be really measured in a general (ideal or 
weak) measurement is the density matrix (with 
all its elements), and only the density matrix! 
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Methods 

An alternative scheme for the tomography of a 
pure state via the weak measurement of a single 
observable. In the main text, we see that it is sufficient 
to perform a weak measurement of a single observable for 
the tomography of a pure-state wave function \ip). This 
single observable could be chosen alternatively as follows. 
We introduce a single observable 

A = ^2\ i \a i )(a i \, (23) 

i 

where the eigenvalues Ai are not degenerated, i.e., A^ ^ 
Xj if i ^ j. We simply perform a weak measurement of 
the observable A and retain all of the data corresponding 
to different final states \bj). When the final state is \bj), 



the weak value of A is given by 

= (bAAty) = E i; (fr> t )A t ( a ^) = EiMi^i 

(24) 

From Eq. (|M|) we have 

- WM^i = 0. (25) 

i 

We define a unitary matrix j3 with the matrix el- 
ements given by the coefficients /3j, and introduce 
the diagonal matrices A = diag{ Xq, Ai, •• • , Xd-i}, 
w = diag{wo,wi, ■ ■ ■ ,Wd-i} and the vector ip = 
(ipOjipi,- ■ ■ ,ipd-i) T , where T denotes the transpose. 
Then, Eq. (|2"5")) can be written as 

(/3X - w[3)~$ = 0. (26) 

With the notation M — j3X — w/3, we have Mnj) = 0. 
Thus, the state tomography in this case is to find the 
vector i/j corresponding to the zero eigenvalue of M. The 
matrix M depends on w, which is determined by the weak 
values that are directly read out from the experiments. 
Due to noise and imperfections in the experiments, ma- 
trix M may not have a zero eigenvalue. Instead, multi- 
plying Mip — by M' from the left, we have 

(M^M)~tj) = 0. (27) 

Therefore, the vector ip is the eigenvector in the kernel 

space of the non-negative matrix M^M. The vector -0 
is uniquely determined if M'M has a one-dimensional 
kernel space. It is not determined, and another set of 
measurements must be performed if M has a kernel 
space of more than one dimension. However, it is also 
possible that the matrix M'M determined from the ex- 
periment may not have a kernel space. In such a case, we 
choose the most likely one, i.e., we choose the vector -0 
as the eigenvector corresponding to the smallest eigen- 
value of M'M. The value of the smallest eigenvalue of 
M indicates the amount of noise and imperfection in 
the experiment. 

Derivation of the position shift and the momen- 
tum shift. Here, we derive the position shift and the 
momentum shift in (f2"0"]) and (f2"Tj) , respectively. The 
time evolution operator corresponding to the interaction 
Hamiltonian in JTSJ) is given by U — e~ l ^< ffi^®P» . If the 
system is successfully post-selected by H,-, which occurs 
with a probability (to the first order of g{) of 

Pj =Tr{{U J ®I)Up s ® p d U^) 

= Tr{U jPs ){l + 2 QilmWji ( Pt )) (28) 

i 

the final (unnormalised) state of the pointers is given by 
(to the first order of g{) 

p' d = Tr s ((n, ® I)U Ps ® Pd Ut) (29) 
= Tr(n ]Ps )( Pd -i^gi(W jiPi p d - W* iPdPi ))(30) 
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Here, Tr s (Tr) denotes the trace over the system (whole) 
Hilbcrt space. The average position shift of the ith 
pointer, conditional on the successful post-selction by IT,-, 

is given by 8qi — T ^f^ d ^ — (<7i). A straightforward cal- 
culation to the first order of gi gives 

Sqi = giReWji + grfmW-ji ({p t - (pi) , q t ~ (©)}) . (31) 



When the initial states of the pointers are well-behaved 
states, such as Gaussian states, the second term on the 
right-hand side of (|3Tj) vanishes, and we immediately have 
(|20p . Similarly the average momentum shift of the ith 
pointer, conditional on the successful post-selection by 
Uj, is given by Sp t = - (Pi), which yields ((21} to 

the first order of gi. 
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